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We report on a study of charmonium at high temperature in 2-flavour QCD. This is the first such 
study with dynamical fermions. Using an improved anisotropic lattice action, spectral functions 
are extracted from correlators in the vector and pseudoscalar channels. No signs of medium- 
induced suppression of the ground states are seen for temperatures up to 1.5r e , while at T ~ 2T C 
there are clear signs of modifications. The current systematic and statistical uncertainties in our 
data, in particular the relatively coarse lattice and small volume, do not allow us to draw a firm 
conclusion at this stage. 
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1. Introduction 

The fate of charmonium states in the deconfined phase of QCD has long been a subject of 
interest, following the suggestion [|I|] that J/y suppression may be a signature of deconfinement 
in heavy-ion collisions. Potential model calculations indicated that J/\j/ might disappear from the 
spectrum almost immediately after deconfinement, providing an unambiguous signature. Recent 
lattice calculations [§, [| |J have cast doubt over this, indicating instead that J/\ff may survive in 
the plasma up to temperatures as high as 1.5 — 2T C . However, all these studies have been earned 
out in the quenched approximation, raising serious questions about their reliability. 

The properties of hadrons in the medium are encoded in the spectral functions p(co,p), which 
are related to the imaginary-time correlator G(T,p) according to 

1 f°° 
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where 
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Here, ns is the Bose-Einstein distribution function. Spectral functions may also be used to extract 
transport coefficients, while at non-zero momentum they contain information about phenomena 
such as Landau damping [ph. 

Determining p(ffl) from lattice correlators G(t) is an ill-posed problem, but it is possible to 
get a handle on it using the Maximum Entropy Method (MEM) It is however crucial to have a 
sufficient number of points in the euclidean time direction, and for this puipose anisotropic lattices 
have a strong advantage over the more common isotropic formulation. A further advantage is that 
it allows the temperature to be varied in small steps while keeping the lattice spacing and spatial 
volume fixed. 



2. Simulation details 

We use the Two-plaquette Symanzik Improved gauge action [^] and the fine-Wilson, coarse- 
Hamber-Wu fermion action ^ with stout-link smearing [g]. The fermion action has been designed 
with heavy quarks in mind, and a quenched study [||] found that the same bare anisotropy can be 
used for valence quark masses ranging from well below strange to well beyond charm. All-to-all 
propagators [[K]] with no eigenvectors and two noise vectors diluted in time, colour and even/odd 



in space, were used to improve the signal. The parameters correspond to point 5 in []1 Using the 
1S-1P splitting in charmonium to set the scale, the lattice spacings were found to be a t = 0.025fm, 
a s = 0.2fm, with the (quark) anisotropy £ ~ 8 determined from the charmonium dispersion relation. 
The sea quark mass corresponds to m n /m p m 0.55. 

Anisotropic lattices with dynamical fermions involve problems that are not present in the 
quenched approximation: since the gauge fields depend on the fermions, it is not possible to first 
fix the gauge anisotropy and then tune the fermion anisotropy so that the physical (measured) 
anisotropy for the fermions matches that of the gauge fields. Instead, the fermion and gauge 



anisotropics must be tuned simultaneously [11]. We do not as yet have data at the fully tuned point 
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given in [11], so unequal quark and gluon anisotropics are still a significant source of systematic 
uncertainty. 

We have performed simulations on TV 3 = 8 3 lattices with N t = 48,32,24 and 16, correspond- 
ing to T pa 0.75,1.1,1.5 and 2.2r c respectively. In all three cases we have used 100 configu- 
rations, sampling configurations every 10 HMC trajectories. In addition we have simulated at 
N t = 28, 30, 33 — 35 in an attempt to locate T c . 

In this preliminary study we only look at zero momentum. The euclidean correlators have 
been computed using local (unsmeared) operators: 

Gr(r) =^{xj/(x,t)rxi/(x,t)xj/(y,t + x)Txi/(y,t + x)), (2.1) 

x,y,t 

with r = 75, Yi, 1,'75^, corresponding to the pseudoscalar, vector, scalar and axial- vector channel 
respectively. The signal in the scalar and axial-vector channels is poor and those results will not be 
shown here. 

The MEM analysis has been performed with the continuum free spectral function ft) 2 as default 
model, using the euclidean correlators in a time window starting at % = 2, and cutting off the energy 
integral in < JTTT] ) at cir&Wx = 6. 

3. Results 

Figure |j] shows the average Polyakov loop (L) as a function of \/N t . We find that with our 
small volume and large lattice spacing it is not possible to determine the critical temperature, since 
no region exhibits a particularly rapid change in this quantity. In agreement with this, the Polyakov 
loop susceptibility has no discernable peak in this region. Larger lattices will be required to deter- 
mine T c . In the absence of any such determination we take T c to be in the region where (L) begins 
to assume a value significantly above zero. 

We have computed correlators for bare 0.2 
quark masses a t mo = 0.1,0.2, where 0.1 is 
close to but slightly lighter than the physi- 
cal charm quark mass. Figure || shows the 
spectral function obtained from MEM for A 
the pseudoscalar (T] c ) and vector (J/ \ff) chan- v 
nels. The position of the main peak for N t = 
48 agrees with the mass obtained for the re- 
spective particles on the same lattices, us- 
ing a variational basis of smeared operators 



0. 1 - 



0.02 0.03 0.04 



[]12[]. The second peak cannot be identified a J 



with the first radially excited state, which on 

these lattices is found to be a t m' = 0.52 [0, Fi 8 ure 1: Avera 8 e Pol y akov loo P as a function of tem " 
in agreement with the PDG value for y'. It P erature - 
is most likely a combination of lattice arte- 
facts and contributions from excited states 2S, 3S etc, which cannot be resolved by the present 
data. 
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The issue of lattice artefacts can be ad- 
dressed by comparing with the free lattice 
spectral function, which is shown in fig. ||, 
together with the continuum free spectral func- 
tions. The most striking feature is the cusp 
at a,(0~ 0.6(0.77) for a t m = 0.1(0.2), which 
coincides with the second peak in our data. 
It is thus not possible to attribute any phys- 
ical significance to the second peak in our 
data. The lattice spectral functions also un- 
dershoot the continuum curve at an early stage, 
indicating that a finer lattice is highly desir- 
able. Due to the lattice cutoff, the free spec- 
tral functions go to zero for a t (0> 1.25, and 
it would thus be sensible to also cut off the 
integral in ( |1.1[ ) near this point. 

Taken at face value, the results shown appear to indicate that the IS states remain unchanged, 
or may even (in the case of r\ c ) become more strongly bound, at temperatures 1-1. 5T C . This would 
not be entirely surprising in light of recent calculations of heavy quark internal energies [|I"3|]. 
However, it will be necessary to investigate to what extent these results remain unchanged with 
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Figure 3: Free-fermion spectral functions in the pseudoscalar and vector channels, for bare quark mass 
a t m = 0.1 (left) and 0.2 (right). 



increased statistics. We find some dependence on the range of time separations used in the analysis; 
this may be statistics-related. 

At N t = 16 the signal has changed substantially: the second peak has disappeared, while the 
primary peak is substantially weakened and shifted to higher energies, indicating a melting of the 
charmonium IS states. There is also a non-zero signal at low CO, but whether this reults in a non- 
zero conductivity requires further investigation. 



4. Discussion and outlook 



The results presented here appear to confirm the picture emerging from quenched simulations, 
that J/y and r\ c survive in the medium up to l.5T c or higher, while excited states begin to melt 
away at lower temperatures. There are even indications that v\ c might be more strongly bound 
at intermediate temperatures, which, if true, would be in accordance with recent results for the 
internal quark-antiquark internal energy [13]. No such changes are seen for J ' j\\r, spoiling this 
agreement. There are, however, a number of issues that must be addressed before any conclusions 
can be drawn with confidence. 

An important source of systematic uncertainty is that the anisotropy has not been completely 
tuned: the quark anisotropy is significantly larger than the gluon anisotropy. This can lead to the 
quarks "feeling" a higher temperature than the gauge fields. It will be important to repeat these 
calculations with a fully tuned parameter set. 

The lattice volume used in this initial study is quite small, at only 1.6 fm 3 , so finite volume 
effects may have a substantial impact on the results. Future simulations will be carried out on a 
larger volume. It is also important to increase the statistics, in particular if we wish to study the 
effects on higher excited states, which will only be clearly resolved with higher statistics. This 
is also necessary if the effects of the zero-temperature decay width (negligible for J /y), thermal 
width and effects of finite statistics are to be resolved. 

As the free spectral functions in fig. || indicate, our coarse spatial lattice means that lattice 
artefacts are substantial even at relatively low energies. Although this problem may be ameliorated 
somewhat by using the free lattice spectral functions as part of the prior knowledge, it will ulti- 
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mately be necessary to repeat the calculation at smaller lattice spacing. This will however require 
a new nonperturbative tuning process. 

These results have been obtained using the continuum free spectral function as default model. 
We are planning to repeat our analysis using other default models, in particular the free lattice 
spectral function, to provide a check on the systematics of the MEM analysis. 

In the future, we plan to extend this study to light vector meson correlators at zero and non-zero 
momentum, which will yield information about dilepton production rates in the plasma. Lattice 
artefacts are expected to be less severe at smaller quark masses, so this may be feasible even on 
current lattices. We also intend to study Landau damping by investigating the behaviour of light- 
and heavy-quark spectral functions at non-zero momentum below the lightcone. 
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